Flavor Oscillation from the Two-Point Function 

Mario Martone 1,2 ' 3, * and Dean J. Robinson 1 '' 

1 Laboratory for Elementary- Particle Physics, Cornell University, Ithaca, N.Y. 
2 Dipartimento di Scienze Fisiche, University of 
Napoli and INFN, Via Cinthia 1-80126, Napoli 
3 Department of Physics, Syracuse University, Syracuse, N. Y. 
(Dated: January 15, 2013) 

Abstract 
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I. INTRODUCTION 



The phenomenon of flavor oscillation plays an important role in the physics of neutral 
meson and neutrino systems. In particular, flavor oscillation provides the only means to mea- 
sure the extremely small mass and decay rate splittings among the neutral mesons, and also 
provides convincing evidence for the existence of non-zero neutrino masses. The theoretical 
descriptions of flavor oscillation fall into several categories, including the basic plane wave 
Pontecorvo formalism [1, 2], intermediate [3-13] and external [14-21] wavepacket approaches 
and quantum field theoretic results [22-31]. Some detailed reviews of these approaches, their 
underlying assumptions, results and difficulties can be found in Refs. [2, 32-34] (and ref- 
erences therein). To be very brief: In the first, one assumes the flavor states are unitary 
combinations of plane-wave mass eigenstates that follow spacetime worldlines, and one must 
carefully define the proper times of the mass eigenstates in order to obtain the well-known 
Pontecorvo oscillation formula. The intermediate wavepacket approach treats the oscillat- 
ing degrees of freedom as a linear combination of one-particle states, while the external 
wavepacket approach treats the oscillating particles as quantum fields, whose propagator is 
convolved with wavepackets at the source and detector. 

A large amount of Literature has been devoted to deriving, studying and comparing 
oscillation formulae within these different approaches. Particular care has been taken to 
include important effects such as measurement uncertainties, coherence effects, the finite size 
of the detector and source, all of which together lead to somewhat complicated formulae. 
Our goal in this paper is less ambitious: Using a quantum field theoretic approach, we 
present a simple formalism of oscillation based entirely on the properties of the spatially 
Fourier transformed propagator. We call this the spatial two-point function. The resulting 
oscillation formulae are particularly elegant, and precisely reproduce both the Pontecorvo 
neutrino result and (CP violating) neutral meson mixing results (see e.g. Ref. [35, 36]) in 
appropriate parameter regimes. 

To construct this formalism, we assume that the oscillation experiment measures the 
exchanged energy E and source-detector displacement L to infinite precision, along with 
flavor at both the source and detector. This is possible because E, L and flavor are com- 
muting observables, so that an amplitude which depends exclusively on these quantities is 
well-defined. The key idea is that the spatial two-point function A(E, L) in the flavor basis 
is a well-defined amplitude which encodes flavor oscillation over a displacement L at energy 
E. We therefore assume that the experiment amplitude is proportional to A(E,L) and 
explore the resulting oscillation formulae. The advantages of this description are: There is 
no ambiguity in the choice of reference frame - all computations are done in the lab frame 
and one never needs to introduce proper times into the formalism; the oscillation probabil- 
ities can be computed exactly; and one obtains formulae whose physical meaning can be 
easily discerned in various limits. Since we do neglect several real physical effects mentioned 
above — in particular the physics of the source and detector are neglected — the limits of 
the applicability of our theoretical description to actual oscillation experiments should be 
carefully examined. Nonetheless, we believe this approach provides an instructive, leading 
order description of the physics of oscillation in real experiments. 

In terms of the previous Literature on this subject, our approach is best categorized as 
a special case of the above-mentioned external wavepacket formalism with stationary states 
[16]. However, to our knowledge, the oscillation physics contained just in the two-point 
function has not been thoroughly investigated and the resulting general oscillation formulae 
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for unstable particles obtained by our approach have not been previously presented. One 
exception is Ref. [22], whose amplitudes for stable fermions agree with our results for the 
special case of stable particles. 

This paper is structured as follows. In Section II we present the oscillation formalism. In 
Section III the exact spatial two-point function A(E, L) for unstable fields is presented, and 
the exact oscillation probabilities and integrated oscillation probabilities are computed. In 
Section IV we examine our results in several different parameter regimes and recover both 
the neutral meson-mixing results and Pontecorvo neutrino oscillation results in appropriate 
limits. 

II. FORMALISM 

A. Experiment Amplitude 

Our starting point is to consider an experiment which involves the propagation between 
a source and a detector of a set of fields {</>"}, which are allowed to mix. As usual, a is an 
experimentally measurable label called the flavor, which is henceforth always denoted by a 
Greek index. The set {4> a } is called the flavor field basis. 

In this paper we assume the a — > j3 oscillation experiment measures the exchanged energy 
E and source-detector displacement L to infinite precision in the lab frame. The amplitude 
for the experiment must then have the form 

M = M a p(E,L). (1) 

This is a well-defined amplitude since E, L and flavor are commuting observables. Note that 
as a consequence of the infinitely precise E and L measurement neither the time of travel 
nor the three-momentum between the source and detector is well-defined, because these 
observables do not commute with E and L respectively. In other words the initial and final 
states of this amplitude must be energy-spatial eigenstates, rather than momentum-time 
eigenstates. 

The key idea of this paper rests on the observation that the time Fourier-transformed 
time-ordered exact two-point function — the spatial two-point function — defined by 

A a p(E,L) = J dt(r{^(t,L)r\0,0)}y Et , (2) 

is the field theoretic object which encodes the oscillation of flavor a — > /3 over a displacement 
L with energy E. (As usual (T , {0 /5 (x)4> a \y)}) = A af3 (x — y) is a function of x — y due to 
translation invariance.) It is therefore natural to write 

M af s(E, L) = A%A^A aP (E, L) , (3) 

(no sum over a, (3) where encode the physics of the source and detector, which we have 
assumed factorizes out of the amplitude 1 . Assuming A$ D are known, the implication of Eq. 



The general criteria under which such a factorization may be possible in real oscillation experiments has 
been examined previously in detail (see e.g. Ref. [34]). Foregoing such a discussion, our intent here is 
that the physics of the source and detector can be neglected up to their ability to distinguish flavor. 
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(3) is that \A a p(E, L)\ 2 is a measurable quantity, from which we may proceed to construct 
oscillation probabilities. 

So far we have not specified the spin of 4> a . As is well-known, if <p a are massive they must 
create spin-j particles, with j a half-integer, which have 2j + 1 spin degrees of freedom. We 
assume that these spin degrees of freedom decouple, so that we need only consider scalar 
propagators henceforth. 



B. Oscillation Probability 

Having written down the amplitude for the experiment, we now define the flavor oscilla- 
tion probability via 

P^(E,L) = l A ^ L )| 2 ( 4 ) 
£|A a7 (£,£)| 2 

7 

Here P a ^p forms a well-defined probability distribution, since P Q( g > and Y2pP a p = 1- In 
some experiments, measurement of |A a/3 | 2 at a precise L is replaced by a volume-averaged 
measurement, 

Ap(E) = J d s L\A a0 (E,L)\ 2 . (5) 

This is equivalent to the time-averaged amplitudes measured in e.g. meson mixing exper- 
iments, in which the initial and final flavor states are determined by tagging via decay 
products (see e.g. Ref. [35, 36]). We can correspondingly define an integrated oscillation 
probability 

^Ug) _ fd?L\A af} (E,L) ~ 
E^a,(E) J2fd3L\A ai (E,L)f 



This is also a well-defined probability distribution. 



C. Propagator and 1PI Basis 



So far in this paper we have formulated a description of flavor oscillation in terms of just 
the exact quantum amplitude A(E, L), which is equivalently defined as the spatial Fourier 
transform of the exact propagator, A{p 2 ). Explicitly, 



&ap(E,L) 



d 3 p 

(2^)3 



A a/3 (p 2 )e 



ip-L 



(7) 



Applying external field methods to the path-integral formulation of quantum field theory, 
it is a well-known result (see e.g. [37, 38]) that for a set of fields {4> a } the exact two- 
point function is the inverse of the exact two-point one-particle- irreducible (1PI) function: 



[x 



y) = 



-n3(* 



y). The Fourier transform of this result is 



A a/3 (p 2 ) = 



p 2 l -M 2 (p 2 ) 



(8) 
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Henceforth we shall call the N x N matrix of functions M 2 (p 2 ) the exact two-point 1PI 
function. 

In general, one cannot compute the exact propagator A(p 2 ) exactly for all p 2 . However, 
the combination of Eqs (7) and (8) suggests that the exact spatial two-point function is 
sensitive only to the pole structure of A(p 2 ). As we shall see below, with suitable assumptions 
this pole structure depends only on physical masses and rest frame decay rates, permitting 
us to construct exact oscillation probabilities in terms of just these measureable quantities, 
E, L, and a mixing matrix, despite our incomplete knowledge of the exact propagator. 

Now, the exact propagator (8) is generally not diagonal in flavor space — there would be 
no oscillation if this were the case — but the analytic structure of A(p 2 ) is greatly simplified 
if the exact propagator can be diagonalized. Ultimately, we want to be able to write 

A a ,( P 2 ) = U^iU-^^ip 2 ) , A,( P 2 ) ee p2 _ l M]{p2) , (9) 

so that 

A a ,(E,L) = U<-(U-y?A 3 (E,L) , A 3 (E,L) ee J ^ (p2) . (10) 

In Eq. (9), U is the constant and possibly unitary matrix that diagonalizes A{p 2 ) (equiv- 
alently M 2 (p 2 )), and M 2 (p 2 ) are the N eigenvalues of M 2 (p 2 ). Below we'll see that the 
Mj{p 2 ) determine the physical masses and rest frame decay rates of the particles propagat- 
ing in A a p(p 2 ). 

In constrast to the usual diagonalization of the classical Lagrangian mass terms, diago- 
nalization of the exact propagator may be non-trivial. In Appendix A we discuss the details 
of the diagonalization of A(p 2 ), the properties of U and how this exact quantum formalism 
both relates to and differs from the usual classical mixing matrix formalism. For our pur- 
poses here, we assume A(p 2 ) is diagonalizable in the manner of Eq. (9). Unless otherwise 
stated, we also assume U is unitary. An immediate consequence of unitarity is that spatial 
two-point function can now be written as 

A afi (E,L) = U ai U f>i *A j (E,L) . (11) 

Let us now define the 1PI basis. This basis is a generalization of the mass basis derived 
in the classical formalism, that may accommodate both unstable particles and a description 
of CP violation for two flavors. In particular, if U is constant (but not necessarily unitary), 
then there exists a well-defined second basis of fields henceforth denoted by a Latin 

index, which are defined via the linear transformations 

Note that $ creates a particle state, while <f> creates an anti-particle state: We have chosen 
this definition of basis change by U in order that it coincides with the usual definition in 
terms of one-particle quantum states. We call the 1PI basis for the following reason. If U 
is unitary, then observe that not only M 2 (p 2 ) but also the two-point function is diagonal, i.e. 
(T{(f> l (x)(f^'(y)}) = 8ijAj(x — y). This implies that M 2 (p 2 ) is the 1PI function for J , whence 
the name. In contrast, if U is not unitary, then even though M 2 (p 2 ) is still diagonalized by 
U, we have (T{0*(x)0- ?t (y)}) ^ 5ijAj(x — y). Therefore M 2 (p 2 ) is no longer the 1PI function 
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for (f) 3 . Nonetheless, we shall always refer to the field basis defined by Eq. (12) to be the 
1PI basis, and often call ft' the 1PI states. 

Tying the diagonalization of the exact propagator and the definition of 1PI basis together, 
we can now explain why we have taken care to consider the case of non-unitary U : We do so 
to accommodate CP-violating two-flavor neutral meson oscillations (for three or more flavors, 
even unitary U may have a CP- violating phase), which we consider in Sec. HID. The idea 
is that the Hamiltonian for such a system is diagonalized by a constant non-unitary matrix 
[35, 36], so we therefore expect U to be non-unitary too. In this context the flavor field 
basis (1PI basis) then corresponds to the CP conjugate states (evolution eigenstates). One 
deduces U is a particular constant 2x2 non-unitary matrix, from which we can immediately 
derive the usual oscillation formulae. 

D. Exact Propagator Analytic Structure 

Let us finally examine the analytic structure of Aj(E,L), which is explicitly 



If ft is unstable, then the propagator Aj(p 2 ) will have a unique Breit-Wigner or resonance 
pole, which by convention is a simple pole located at 



Here rrij is the physical mass and Tj > is the rest frame decay rate. The non-zero 
imaginary part for this pole enforces the usual Feynman pole prescription and associated 
time-ordering, so that we need not add the usual ie convergence term in the denominator 
of Eq. (13), provided we assume Tj ^ 0. Consequently, taking the Tj — > + limit, which 
corresponds to ft being stable, can only be performed after all integrations and other limits 
are evaluated. 

By definition there are no higher order poles in Aj(p 2 ) and the residue at the pole (14) is 
unity: Eq. (14) and this latter condition are equivalent to M 2 {m 2 — irrijTj) = m 2 — irrijTj 
and M 2 '{m 2 — irrijTj) = respectively. 

III. EXACT OSCILLATION PROBABILITY 
A. Spatial Two-Point Function 

Computation of the oscillation probabilities (4) and (6) boils down to computing the 
spatial two-point function Aj(E,L). As shown in Appendix B, the integral (13) can be 
performed exactly, with final result (BIO) 




(13) 



p 2 = m 2 — irrijTj . 



(14) 



Aj(E,L) 



i 



[ % r / i 1/2 i 



R) + A 2 - Rj L 



nl/2 



(15) 



in which 



Rj = E 2 — m) 
Aj = rrijTj . 



■j ' 



(16) 
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Note that the exact result in Eq. (15) is independent of the orientation of L 



B. Exact Probabilities 

We may now compute the exact oscillation probability via application of Eqs. (4), (10) 
and (15), and the exact integrated oscillation probability via Eqs. (6), (10) and (15). It is 
convenient to define the wavenumber and characteristic inverse decay lengths 



Uj = 



1 

i 

72 



y/H* + Aj + Rj 



,1/2 



R] + A]- Rj 



1/2 



along with 



Auj k = uij — u k , AQ k = Q — ( k , ( jk = (j + (, 



k ■ 



We call Auj k the oscillation wavenumber. 

Exploiting the unitarity of U, one finds the exact oscillation probability 



P a ^(E,L) = <j J2\ Uaj \ 2 W Pj \ 2 e^ L + 2j2~Re[u aj U l3j *U ak *U l3k t 

j<k 

-l 

j 



iAujkL — CjkL 



^2\U aj \ 2 e~ 2CjL 



(17) 



(18) 



(19) 



For Q — > + , this has the exact form of the Pontecorvo oscillation formula [1, 32]. We will 
show below that within a certain parameter regime, Auj k ~ (ml — m 2 )/2E, recovering the 
usual result. The exact integrated oscillation probability is similarly 



j<k 



jjajjj/3j*jjak*jjl3k Cjk + iAuJfh 



X 



J2\U aj \ 2 /2Q 



-i -l 



(20) 



(The Aj normalization i/iiiL plays an important role in computing the integrals in Eq. (6).) 
Note that for both Eqs. (19) and (20) we have not assumed CP conservation. 



C. Two-Flavor Formulae 



It is particularly illuminating to present the oscillation probability and integrated oscil- 
lation probability for the case that there are just two flavors. In this case we can choose U 
to be real, orthogonal: The only physical parameter is the mixing angle and there is no CP 
violation. We adopt the convention for two flavors that a = +, — and j = 1,2. One obtains 
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oscillation probability 



P a ^(E,L) = ! [ \U al \ 2 \Ue i \ 2 
+ {\U a2 \ 2 \U^ 2 \ 2 



\U al \ 2 + \u a2 \ 2 e~ 2A< * 2lL 



+ <^ 2U al U l31 U a2 U^ 2 cos(Aw 12 L) 



\jjali2 e +AC 2 iL + ijja2i2 e -A( 2 iL 



• (21) 



Assuming without loss of generality that A( 2 i > 0, then for A£ 21 L ^> 1, the first term is 
asymptotically constant, the second decays to zero while the third term produces a damped 
oscillation decaying to zero, with oscillation wavenumber Ao>i 2 . In particular, for A( 2 iL ^> 1, 



P^{E,L)~\U^ 



(22) 



If we adopt the notation that (0^)^ creates \f3) and (0 J )^ creates \j), then the right side of 
Eq. (22) is nothing but |(/3|1)| 2 . This is the probability of the (3 flavor state being measured 
as the j = 1 1PI state, which has the longer decay length 1/d > l/( 2 . This behavior is 
familiar to that found in the K neutral meson system: Since the Kl eigenstate has a much 
longer decay length than the K$, then the Kl will be exponentially more abundant at large 
distances from the source compared to the K$ state. As a result, at large distances there 
is no more oscillation and the oscillation probabilities K — > K or K — > K both collapse to 
\{K\K L )\ 2 . This is exactly the behavior in Eq. (22). 

Before presenting the two-flavor integrated oscillation probability, for convenience we first 
define 

A ^12 _ A( 2 i . . 

V = ~T— ■ (23) 



x 



12 



c 



12 



In Sec. IV B below we shall verify that x and y reduce to their usual definitions x ~ Am/f 
and y ~ Ar/2f within a certain regime of the parameters E, m 12 and r 12 . With the 
definitions (23), Eq. (20) reduces to 



P l a ^(E) = \ Ifr 1 ! 2 !^ 1 ! 2 ^ + y) + |f/ a2 | 2 |f^ 2 | 2 (l -y) + 2U al U^ 1 U a2 U (32 
l+y(\U al \ 2 -\U a2 \ 2 ) 



i-y 2 

l + x 2 



X 



(24) 



D. CP Violation and Non-Unitary Diagonalization for Two Flavors 

The above oscillation probabilities (19) and (20) (or (4) and (6)) may be generalized to 
the case that U is constant and non-unitary, which is applicable to the study of CP violation 
in two-flavor neutral meson mixing. In this context, we identify the flavor fields (0 ± )^ as 
the creation operators of the CP conjugate states \P°) and \P°), while the 1PI basis ((j) 1 ' 2 )^ 
create the evolution eigenstates \Pl,h) respectively. Comparing Eq. (12) with the usual 
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notation for \P°) and \P°) in terms of \Pl,h) (assuming CPT symmetry) 

\P°) = ^(\Pl) + \Ph)), 

\P*) = ^{\Pl)-\Ph)) , (25) 

then leads to the identification 



U = — + ( q q ) , (26) 
2pq -\p -pj y ' 

which is non-unitary for \p/q\ 7^ 1. As \p/q\ 7^ 1 is sufficient for CP violation in Eq. (25), the 
consequence of the identification (26) is that CP violation in two-flavor mixing is manifested 
as non-unitary diagonalization of the two-point function. Note also that if U is non-unitary 
then the 1PI states are no longer orthogonal, as expected for the evolution eigenstates \Pl,h)- 
That is (Tcf) l (x)(f>>^(y)) 7^ SijAj(x — y). From Eq. (9) we have for non-unitary U 

A a/3 (E,L) = U aj (U~ 1 ) jl3 A j (E, L) , (27) 

and from Eqs. (15) and (17) one then finds exact amplitudes 

|A±±(£, L)\ 2 = —l-^e-^ 1 [cosh(ACi2^) + cos(Au; 12 L)] , 
\A + _(E,L)\ 2 = ^_l- i e-^ i [cosh(AC 12 L) - cos(Ao; 12 L)] 



A_ + (E,L)\\ (28) 



These strongly resemble the amplitudes found from the usual meson mixing quantum me- 
chanical analysis (see below), except that here a spatial dependence has replaced the usual 
time dependence and Aujk, AQk an d (jk have the form presented in Eq. (18). One finds 
exact integrated oscillation probabilities 



p i _ 2 + x 2 - y 2 



2 + x 2 — y 2 + \q/p\ 2 {x 2 + y 2 ) 



PU-= : „ ," 2 ^ — - (29) 



x 2 + y 2 + \p/q\ (2 + x 2 — y 1 



and one can also contemplate measuring the ratio of the amplitudes for oscillation into either 
flavor state 

A\_(E) _ fdL\A + _(E,L)\ 2 
A\ + (E) JdL\A ++ (E,L)\ 2 
2 2 1 2 

= 1 x + y m 

p 2 + x 2 - y 2 ' V } 
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Let us compare the exact results in Eqs. (29) and (30) with the analogous formulae ob- 
tained via the usual quantum mechanical treatment of neutral meson oscillations. Following 
Refs. [35, 36] we can write down the time evolution for an initial pure \P°) and \P°) state 



\P a (t))=g + (t)\P ) + ±g-(t)\I*>) 

V 

mt)} = ^9-(t)\p o }+ 9+ (t)\n , 



where 



,-rt 



cosh(Art/2) ± cos(Amt) 



In the standard notation 





x = Am/r , 


y 


= AT/2T , 


where T = (Ti + T 2 )/2, the two formulae in Eqs. 


(29) should be compared with 




f~dt \(P°\P°(t))\ 2 




o i 2 2 

2 + x -y 


ST** 


\(P°\P°(t))\ 2 + \(PV\P°{t))\ 2 




2 + x 2 -y 2 + \q/p\ 2 (x 2 + y 2 ) 


and 






x 2 + y 2 




(p°\p°m 2 + \(p o \p°(t))\ 2 




x 2 + y 2 + \p/q\ 2 {2 + x 2 -y 2 ) 



(31) 

(32) 
(33) 

(34) 
(35) 



Finally Eq. (30) should be compared to 

f °° dt \(p°\p°(t))\ 2 



x 2 + y 2 
2 + x 2 -i 



(36) 



Our exact results are in perfect agreement with those of the usual analysis, except that in 
Eqs. (29) and (30), the parameters x and y have the more general definitions encoded in 
Eqs. (16), (17) and (23). Once again, in Sec. IVB below we shall verify that x and y 
reduce to their usual definitions x ~ Am/f and y ~ AT/2T within a certain regime of the 
parameters E, and r^. 

Before proceeding, please note that the oscillation probabilities (19), (20), (21), (24) and 
(28) presented in this section are a function of only Aujk, AQk, (jk or Q (of the latter three 
variables, only two are independent). Consequently, specifying just u>j and Q is sufficient to 
specify the oscillation probabilities. This shall be our practice throughout the remainder of 
this paper. 



IV. REGIMES 

The results presented in Sec. Ill are elegant and concise, but their physical interpretation 
is not obvious. However, in various regimes of the parameters E, rrij and Tj, our exact 
results for the wavenumber and characteristic inverse decay lengths u and ( reduce to simpler 
expressions with clear physical meanings. In this section we explore several different regimes 
of physical interest, and show that in certain regimes our results reproduce the well-known 
neutrino and neutral meson oscillation formulae. 
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A. Particle Regime 



The first regime of interest is the case 

Rj > Aj , i.e. E 2 - > rrijTj . (37) 

It is straightforward to expand Eq. (17) about Aj = 0, and to leading order in Aj/Rj one 
obtains the wavenumber and characteristic inverse oscillation lengths 

2^E 2 -m 2 



to, ~ JE 2 - m) , Q - , [ 1 „ ■ (38) 



The oscillation probabilities (19) and (20) follow immediately from this and Eqs. (18), as 
do their two-flavor versions (21) and (24). In particular, in this regime we have to leading 
order in Aj/Rj 

A 1 /^R~ l + A 2 /^R~ 2 ' A l /y/R~ l + A 2 /^R~ 2 ' 

In this regime, the spatial two-point function for a 1PI state 



Aj(E,L) ~ — - expliJE 2 - m 2 L . 3 3 L > . (40) 

47rL I V 2^E 2 -m) J 

For Tj — > + , this looks precisely like the propagator of an on-shell one-particle state with 
momentum p = (E 2 — m 2 ) 1 / 2 . We therefore call the regime (37) the particle regime. The 
resemblance of the amplitude (40) to that of a particle suggests that we should obtain both 
the neutrino and meson mixing oscillation formula within the particle regime. We explicitly 
verify this in the next two sections. 

The physical meaning of the spatial two-point function perhaps becomes more clear if we 
define analogous Lorentz factors and proper time 



-y^E/mj, Pj = yjl- m 2 /E 2 , Tj = L/^fr . (41) 

Substituting these into Eq. (40) we obtain the spatial two-point function in terms of Tj 
instead of L, which we can interpret as a 'rest frame' propagator. Explicitly, 

Aj{Tj) - i^i exp { imj{ ^ " 1)Tj " • (42) 

The second term in the exponential looks like the usual rest frame decay of an unstable 
particle, and in particular it is clear that Tj can be interpreted as the rest frame decay rate. 
The first term looks like the usual proper time evolution of a particle, except for the 7 2 — 1 
factor. This factor arises because pL is not a Lorentz invariant, but rather E^t — pL = mr 
is. It is a consequence of the experiment measure E rather than the time of transit between 
the source and detector. 

Let us now proceed to verify that the usual neutrino and neutral meson mixing oscillation 
formulae are obtained in this particle regime. 



11 



B. Small Mass Splitting: Neutral Meson Oscillation 



In all known neutral meson systems, the mass difference between the two mass eigenstates 
is extremely small in comparison with their masses. For the K, D, B d and B s neutral meson 
systems one finds [32] 



/ Am 
\ m 



K 



Am 



m 



10 



-14 



D 



[Am\ 
V m J 



10 



-13 



B d 



Am\ 
m J 



10 



-12 



(43) 



It seems then, that the appropriate regime for neutral meson oscillation is the particle regime 
with the additional constraint that the mass splitting is small. We define the mean mass 
m and mass splitting Am via mi = m + Am/2 and m 2 = m — Am/2, so the small mass 
splitting limit is Am/m <C I. We also define y = Ar/2r, in which f = (Fi + F 2 )/2 and 

Ar = r 2 -rx. 

Expanding the particle regime expressions (39) for x and y in the small mass splitting 
limit is complicated by the fact that neither x nor y can be expressed as function of Am/m 
alone. However, one may show that 



x ~ 



Am 

Ar 

2F 



1 + 



1 + 



yo_Am _ 

2(3 2 m 

1 — 2/q Am 

2(3 2 y m 



^ X p (y ,m,E) f Am \ F 

p—2 



p=2 



Y p (y ,m,E) / Am 
2p/3 2 p V m 



(44) 



for which X p and Y p are rational functions of y , m and E. The parameter j3 is defined 
as in Eqs. (41), but for mass m. In general, any of X p , Y p or l/j3 could be arbitrarily 
large for some configuration of the parameters yo, m and E, so the expansions (44) are not 
always well-controlled power series in Am/m. However, it's plausible that the expansions 
are well-controlled in the parameter space regimes, and respective very small mass splittings 
(43), relevant to the neutral meson systems. In such regimes, we then obtain for sufficiently 
small mass splittings 

Am ^ ,AZ\ 

x - — > y - 2f • ( 45 ) 

These are precisely the usual definitions for the parameters x and y in the neutral meson 
mixing formalism. Moreover, in terms of the 'proper time' r - defined for m in Eqs. (41) - 
the same expansion in small mass splitting renders the CP violating amplitudes (28) 



\A ++ (E,L)\ 2 ^ 
|A + _(£,L)| 2 ^ 



32tt 2 L 2 



exp(— Ft) cosh(Arr/2) + cos(Amr) 



32tt 2 L 2 H 2 



exp(— Ft) cosh(Arr/2) — cos(Amr) 



(46) 



Up to a normalization of l/lQir 2 L 2 , these are exactly the time evolution amplitudes found 
within the usual meson mixing analysis [35, 36] as is evident in Eq. (32). The physical 
interpretation of r is the proper time elapsed in the rest frame of a classical particle with 
mass m and lab frame energy E that traverses a distance L. 

Since we have already verified that, in terms of x and y, our integrated oscillation proba- 
bilities match those found in the usual treatment, we have therefore recovered the well-known 
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meson mixing amplitudes and time-integrated probabilities from the structure of the two- 
point function alone. Our analysis, however, also implies that the usual meson mixing results 
are valid only within the small mass splitting particle regime. Outside this regime the more 
general results of Sec. Ill D will apply. 

An immediate question is whether the regimes of validity of our derivation and the stan- 
dard quantum mechanical one disagree. The standard derivation is performed in terms of 
time evolution, and requires a common proper time for the mass eigenstates [35, 36], which 
are one-particle states. If the energy of the oscillation experiment is fixed, as we assume 
throughout this paper, then this assumption is equivalent to assuming Am — > 0. So the 
standard derivation is similarly applicable only in the small mass splitting limit. 



C. (Ultrarelativistic) Stable Particle Regime: Neutrino Oscillation 

For a neutrino oscillation experiment we expect the neutrinos to be ultrarelativistic and 
stable in the lab frame. The ultrarelativistic stable limit of the particle regime corresponds to 
E ^> rrij and Tj — > + for all j. Expanding in this limit, the wavenumbers and characteristic 
inverse decay lengths (38) become, to leading order in m,j/E, 

o 

Tfl ■ 

u^E-^L, = 0, (47) 



so that 



A J (E,L)^—expUEL- i ^LL}. (48) 



Applying Eqs. (19) and the unitarity of U leads immediately to 

T f Am 2 

P a ^(E, L)~J2 \U aj \ 2 \UP j \ 2 + 2 Re U aj U^ j *U ak *U^ k exp I - i^p^L 

j j<k L ^ 

r i /Am 2 \ 

= 5 a(S + 2j2 lm UJ a] U^*U ak *U pk sin I —L J 

j<k ^ ' 

- 4 Re [u aj U f}j *U ak *U f}k '\ sin 2 (^rL) , (49) 
j<k ^ ' 

where Am 2 - k = m 2 — m\. This is precisely the Pontecorvo neutrino oscillation formula. 
Hence we have derived the neutrino oscillation formula in a purely quantum field theoretic 
formalism, involving just the structure of the spatial two-point function. Comparing Eqs. 
(38) and (47), it is straightforward to generalize this result to just the stable particle regime 
E > rrij, Tj — > + via the replacement in Eq. (49) 

/\ 777, ^ I j 

~* ' ' -)• sJE 2 - m 2 - yjE 2 - m\ . (50) 



2E 



D. (Deep) Virtual Regime 



Having verified that our exacts results reduce to the expected results for both the neutral 
meson and neutrino systems, let us now exploit the generality of Eqs. (15), (17) and (19) 
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to push E, rrij and Tj into non-standard, though physically relevant, regimes of parameter 
space. So far we have only considered regimes for which E > rrij, so let us now consider the 

case 

-Rj > Aj , i.e. m 2 - E 2 > rrijYj , (51) 

which we call the virtual regime for reasons outlined below. If also E <C rrij, then we call 
this the deep virtual regime. As we will investigate in Sec. IV E below, the virtual regime 
is particularly interesting if one 1PI state is very heavy compared to another. 

In the virtual regime, the wavenumber and characteristic inverse decay lengths become 
to leading order in Aj/\Rj\ 



^ P^— , Cj^Jrrvj-E*. (52) 
2^m 2 - E 2 V 

Again, the oscillation probabilities (19) and (20) follow immediately from this and Eqs. 
(18). Within this regime, two flavors with a sufficiently small mass splitting have integrated 
oscillation probability described by the parameters 

Ar Am 

x - ' v - ~y~ ■ ( ) 

These are, of course, just a swap of the usual parameters one sees in the particle regime. 
In the virtual regime Eq. (15) becomes 



Aj(E, L) ~ —— exp I i ™ jTj L - Jm 2 - E 2 L } . (54) 
^ L I 2^t~^E 2 V 3 

The spatial two-point function no longer looks like that of a one-particle state. This is 
especially clear in the stable virtual case rrij > E and Tj — > + , for which Aj(E,L) is 
just an exponential decay: This was noticed previously in Ref. [22]. Note also that for the 
unstable case, the wavenumber is determined by the decay rate, rather than by a momentum 
(E 2 — m 2 ) 1 / 2 , and vice versa for the characteristic inverse decay length. 

Let us briefly comment on the physical interpretations of the virtual regime results, from 
which we derive its name. As explained in Sec. II A, the spatial two-point function Aj(E, L) 
does not encode the propagation of just a single particle with a definite momentum. Rather, 
as suggested by Eq. (13), we may think of the spatial two-point function as the contin- 
uous sum of a set of propagators, each corresponding to the propagation of a momentum 
eigenstate. The condition E < rrij then implies all these momentum eigenstates must be 
off-shell, so in this case Aj(E, L) includes no on-shell propagating particles. That is they are 
virtual particles, whence the regime name. Alternatively, E < rrij is analogous to the usual 
quantum mechanical tunnelling condition, with the mass acting as the potential barrier. 
From either point of view, we emphasize that we should expect A(E, L) to be exponentially 
suppressed in the stable case, precisely as we see in Eq. (54). 



E. Mixed Regime 

It is interesting to consider the case that different 1PI states occupy different regimes. 
For example, we could consider a two-flavor oscillation in the case that one 1PI state is in 
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FIG. 1. Quark oscillation experiment t — > Xp. 



the particle regime, while the other is in the virtual regime. We call such a case the mixed 
regime. This scenario doesn't occur for the neutral meson or neutrino systems because the 
mass splitting between 1PI states is very small compared to E. However, the large mass 
hierarchy of the quark sector combined with the possibility of quark oscillations [39, 40], 
provides a natural setting in which we may contemplate a mixed regime oscillation. 

For concreteness, let us suppose there is a fourth quark doublet (t',b'), with masses 
much larger than the top quark mass. The existence of a fourth quark family is strongly 
constrained by the electroweak weak precision measurements [41-43], but nonetheless phe- 
nomenologically still perfectly viable (see Ref. [32] for mass bounds). A quark oscillation 
experiment could then involve a top quark decaying to a final state t — > Xp via intermediate 
b or b' down-type quarks. The generic diagrammatic form of such an experiment is shown 
in Fig. 1. 

Let us adopt the following notation. The flavor of the down-type quarks is determined 
by their up-type partner, so we denote the down-type flavor quarks by b t and b t >. That is, 
a = t,t'. Correspondingly the 1PI states are denoted b and b', so j = b,b'. The idea here 
is that the top quark produces the flavor quark b t at the source vertex, S, while the generic 
final state Xp in the detector can tag the flavor at vertex D. In order to describe the physics 
of this experiment using our formalism, and for simplicity, we also assume the following: 

i) The amplitude of the experiment is described by Eq. (3). 

ii) We neglect the presence of the other two down-type quarks d and s, and consider an 
effective two-flavor mixing between the third and fourth quark generations. Conse- 
quently, the final state Xp only measures flavors (5 = t, t' . 

iii) The b is in the stable particle regime (E > nib and Tb — > + ), while b' in the stable 
virtual regime (E < my and Ty —> + ). 

iv) The energy, E, exchanged between S and D can be precisely measured. 

v) The 2x2 mixing matrix U, which diagonalizes the 1PI function, is unitary. 

The extent to which these assumptions are applicable to an actual quark oscillation experi- 
ment is questionable. Our intent is merely to demonstrate that with such assumptions, we 
can perhaps gain insight into the physics of quark oscillations by use of our formalism. 
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With these assumptions, we have wavenumber and characteristic inverse decay lengths 



u b = sjE 2 -m 2 , u v = , ( b = , Cv = yj rr% - E 2 , (55) 

so that the oscillation wavenumbers 

Au w = \jE 2 -m\ , A( b , b = \Jm 2 -E 2 . (56) 

The oscillation probabilities follow immediately from Eqs. (21) and (24), while the corre- 
sponding spatial two-point functions 

A b (E, L) = -^ L exp ^E 2 -m 2 b L^ , 

&v(E,L) = ^exp | - yJml-E^ . (57) 

In particular, note that the b' two-point function is exponentially suppressed, as we expect 
for a virtual particle. Further, the integrated probability has parameters 



l E 2 


-m 2 b 


ml 


— E 2 



so that in the two-flavor integrated oscillation probability (24) there is no longer any inter- 
ference term — and hence no oscillation — between the two mass eigenstates. From Eq. 
(20) one finds integrated oscillation probability 

r t ^\E) - 2 1+ \ ( jtb\2_\ U tb>\2 ■ ^> 

(Note that convergence of the integrated amplitude with = is ensured by the usual 
ie term, which is taken to zero after integration. Equivalently, since mbTb acts as the e 
throughout this paper, the stable limit is determined by taking Tb — > + after integration 
over L.) As expected, the probability is controlled purely by the mixing of the flavor states 
with the particle-like 1PI state b which is in the stable particle regime. 



F. Threshold Regime 

One last regime of interest, which to the knowledge of the authors has not been previously 
discussed, is the case 

|.Rj| < Aj , i.e. \E 2 - m]\ < mjTj . (60) 

We call this the threshold regime, since E ~ mj. To zeroth order in \Rj\/Aj, the wavenumber 
and characteristic inverse lengths reduce to 

^-Q-^S- (61) 
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and the oscillation probabilities follow as usual. This time (15) becomes 



A,(£,L)^exp|(.-l)^&^ . (62) 

Here, curiously, the inverse decay length and wavenumber both depend on the geometric 
mean of the decay rate and mass, and coincide. We are unaware of an intuitive physical 
reason why they should coincide at threshold. There is, however, a limited particle analog 
to this behavior. If we were to interpret uoj as the momentum, as we did in the particle 
regime, then we would have 

p)^m)-m 3 T 3 /2. (63) 

That is, the 1PI states can be thought of as virtual particles slightly perturbed from the 
mass shell if Tj <C rrij. 

In the threshold regime, a small mass splitting for two flavors results in 

y/T\ — y/T^ 

x = y~ (64) 

V 1 1 + V 1 2 

while if also the decay rates have a small splitting, 1^2 = T =p Ar, Ar <C T, then 

i = ?/~-Am/m-Ar/r«l. (65) 

A well-motivated example of oscillation in which the threshold regime is applicable to both 
1PI states is unknown to the authors. Despite this, we do wish to emphasize that the gen- 
erality of Eqs. (19) and (20) permits exploration of parameter regimes in which a quantum 
mechanical treatment might be unfeasible. 



V. CONCLUSIONS 

In this paper we have used only the structure of the spatial two-point function A a p(E, L) 
to derive general flavor oscillation probability formulae for unstable fields. We have not 
only shown that this structure reproduces the usual Pontecorvo neutrino oscillation formu- 
lae and time-integrated (CP violating) neutral meson mixing formulae, but we have also 
found generalized exact expressions with natural physical interpretations in several different 
parameter regimes. Our results for the stable particle and stable virtual regimes agree with 
the results of Ref. [22] for stable fermions. However, our exact oscillation probabilities for 
unstable fields and the analysis of the unstable particle, threshold and virtual regimes has 
not been previously presented. 

The advantages of the formalism we have employed in this paper are several. The exact 
computability and integrability of A a p(E,L) permitted us to obtain exact, elegant proba- 
bility oscillation formulae. Moreover, the choice of reference frame throughout this paper 
is the unambiguous laboratory frame: There is no need in our approach to contemplate 
mass eigenstate rest frames and proper times. To the extent that complicating effects such 
as coherence, finite detector and source size, non-trivial source and detector physics, and 
measurement uncertainty can be neglected, our results provide an instructive leading order 
description of the physics of flavor oscillation, that is valid over the entire E, m, T parameter 
space. 

In terms of future work, keeping in mind the large existing Literature on this subject, 
perhaps the most interesting avenue left to explore is the analogous formalism for flavor 
oscillation in matter, that is, the Mikheyev-Smirnov-Wolfenstein effect. 
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Appendix A: Diagonalization of the exact propagator 

In this appendix we discuss the subtleties involved in diagonalizing the exact propagator 
in Eq. (8), 

(Al) 

In general the exact propagator A(p 2 ) is not necessarily Hermitian. It is therefore not 
always diagonalizable by a unitary matrix and may not even be diagonalizable at all. If, 
however, A(p 2 ) is diagonalizable by some invertible matrix U, then observe that: We should 
generally expect U to be a function of p 2 , U = U(p 2 ), since A = A(p 2 ); The diagonalizability 
of A a/3 (p 2 ) is equivalent to that of M 2 g(p 2 ), since if one is diagonalizable by U{p 2 ) then so 
is the other. 

Keeping these two observations in mind, in the standard field theoretic oscillation formal- 
ism one first diagonalizes the tree-level Lagrangian mass terms, thus obtaining free propa- 
gators for the mass eigenstates, and then one can construct two-point amplitudes perturba- 
tively. For example, the well-known PMNS (CKM) matrix diagonalizes the lepton (quark) 
masses in the Standard Model with right-handed neutrinos (SM + vr). However, in general 
such a diagonalization does not persist to all orders in perturbation theory. In particular, 
in the SM + vr model the flavor changing 1PI functions are zero at tree-level, but receive 
non-zero contributions at loop level, which are small due to the GIM mechanism. To see 
this, note that the exact propagator for a left-handed neutrino in the mass basis is 




in which we use the usual SM notation. The mass splittings of the leptons ensure that the 
neutrino exact propagator is not diagonal at all loop orders in the PMNS basis. Another 
manifestation of this effect is that flavor changing neutral currents do not appear at tree 
level, but they do appear at higher loop orders. The moral is that if the propagator is 
diagonalizable in the exact theory, it is generally diagonalizable by a p 2 -dependent matrix, 
which is different from the PMNS or CKM matrix at subleading order and not necessarily 
unitary. 

Due to the GIM suppression, it is common in oscillation formalisms to neglect this effect 
because it occurs at subleading order in perturbation theory. Instead, one presumes that the 
propagator is diagonalized by the constant, unitary PMNS or CKM matrix. In this paper 
we make a similar assumption in Eq. (9) in the main text. The validity of this assumption 
is model dependent, and a discussion of the general circumstances under which it applies is 
beyond the scope of this paper. Nonetheless, as the above SM + Vr example demonstrates, 
it is true at leading order in perturbation theory for certain important theories. 



A a p{p 2 



p 2 l-M 2 (p 2 ) 
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We emphasize finally that rather than arising from a diagonalization of the bare classical 
Lagrangian, the matrix U here diagonalizes the exact propagator A a/ 3(p 2 ), which includes 
all quantum corrections. For the SM + vr example discussed above, U therefore coincides 
at zeroth order with the neutrino PMNS or quark CKM matrix. However, in general U acts 
as the mixing matrix between the flavor field basis and the 1PI basis (defined in the main 
text), rather than between the flavor basis and the mass basis of the classical Lagrangian. 



Appendix B: Computation of Aj(E,L) 

In this appendix we compute the integral in Eq. (13): 



^•WlMw (B1) 

First, it is convenient to partition the 3- momentum as 

p = p L L/L +p ± , p ± -L = 0, (B2) 

so that p - L = p L L and d 3 p = d 2 p±dp L . It is clear that the integrand of Eq. (Bl) has a p L 
pole determined by Eq. (14), which becomes in terms of pl 

p 2 L = E 2 - p\ - m 2 + irrijTj . (B3) 

Since L > 0, one can close the pi integration contour on the upper-half complex plane, and 
then only the pi pole in the upper-half plane contributes to the integral. Let this (positively 
oriented) integration contour be denoted by C. (Note that a M 2 (p 2 ) branch cut on the pi 
real axis doesn't affect the integral, since it can be rotated off the axis by an appropriate 
choice of the principal branch.) With the notation of Eq. (B2) the spatial two-point function 
becomes 

"d 2 p± j dp L ie iphL 



L) = ]W? f c ^E 2 -p 2 ± -p 2 L -m 2 + i m 3 T 3 " (B4) 

Note that Aj(E,L) is independent of the orientation of L, so Aj(E,L) = Aj(E,L). The 
physical consequence of the pl contour integration is to force the 4-momentum of the in- 
tegrand (13) to be on the 'pole shell' in the complex sense defined by Eq. (14). We can 
interpret the remaining d 2 pj_ integral to be a sum over on-pole-shell transverse 3-momenta. 

Let us now perform the Pl integral. In more compact notation, the propagator has a. Pl 
pole satisfying 

pl(p±) = R j (p±)+iA j , (B5) 

with 

R j (p 1 _) = E 2 -pl-m 2 , 

Aj = m l V l . (B6) 

We do not make any assumption regarding the sign of Rj(p±). However, for Tj > it is 
clear that Arg[p|] G (0, 7r]. Therefore, defining z = Aj/Rj(p±), it must be that 

Arg[p|] = Tan- 1 ^) ee ^ z > ° , (B7) 

7r — tan \z\, z < 
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where tan -1 | • | : [0, oo) — > [0,7r/2]. This permits us to compactly write p 2 L {p±) in complex 
polar notation. Taking a square root is now trivial, and the pi pole in the upper-half complex 
plane is 



Pl{p\ 



R 2 Ap±) + A 



1/4 



cxp 



^Tan- 



(B8) 



Applying the residue theorem to Eq. (B4), and observing that Pl(p±) is only a function of 
the magnitude of p±, we now have 



A j (E,L) = 



2 p ± dp 



- 1 £ iPL(p±)L 



4tt7 Pl(p±) 

Observe, furthermore, that Eq. (B5) implies p±/pl(p±) = — Pl(p±)- Hence 



4tt Jo 



dp ± p' L (p L )e ip ^ L 



1 Jpl(Q)l 
AuL 



as Pl(p±) — >■ 200 when pj_ — > oo. Writing Rj = Rj(0) = E 2 — m|, we have finally 

1/4 



A 3 (E,L) 



i 



exp <^ / 
exp 



cxp 



2 m 1 Mi 
I™ 1 U 



L 



72 



R 2 + A 2 + Ri 



1/2 i 

L ~T2 



R 2 + A 2 - IL 



1/2 



(B9) 



(BIO) 



In the last line we have used several trigonometric identities along with the definition of 
Tan' 1 in Eq. (B7). 
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